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We derive the Landau-Lifshitz-Bloch (LLB) equation for a two-component magnetic system valid 
up to the Curie temperature. As an example, we consider disordered GdFeCo ferrimagnet where 
the ultrafast optically induced magnetization switching under the action of heat alone has been 
recently reported. The two-component LLB equation contains the longitudinal relaxation terms 
responding to the exchange fields from the proper and the neighboring sublattices. We show that 
the sign of the longitudinal relaxation rate at high temperatures can change depending on the 
dynamical magnetization value and a dynamical polarisation of one material by another can occur. 
We discuss the differences between the LLB and the Baryakhtar equation, recently used to explain 
the ultrafast switching in ferrimagnets. The two-component LLB equation forms basis for the large- 
scale micromagnetic modeling of nanostructures at high temperatures and ultrashort timescales. 

PACS numbers: 75.78.Jp, 75.40. Mg, 75.40.Gb 



I. INTRODUCTION 

The Landau-Lifshitz-Bloch (LLB) dynamical equation 
of motion for macroscopic magnetization vector [1] has 
recommended itself as a valid micromagnetic approach at 
elevated temperatures [2], especially useful for tempera- 
tures T close to the Curie temperature Tc (T > 3Tc/4) 
and ultrafast timescales. In several exciting novel mag- 
netic phenomena this approach has been shown to be a 
necessary tool. These phenomena include laser-induced 
ultrafast demagnetization [3]-[6] , thermally driven domain 
wall motion via the spin-Seebeck effect [7], spin-torque ef- 
fect at elevated temperatures [HIE] or heat-assisted mag- 
netic recording [TU] , 

In the area of laser-induced ultrafast demagnetization, 
the LLB equation has been shown to describe adequately 
the dynamics in Ni [5] and Gd [6]. The main feature of 
the LLB equation allowing its suitability for the ultra- 
fast magnetization dynamics is the presence of longitu- 
dinal relaxation term coming from the strong exchange 
interaction between atomic spins. Because the exchange 
fields are large (10—100 T), the corresponding character- 
istic longitudinal relaxation timescale is of the order of 
10-100 femtoseconds and thus manifests itself in the ul- 
trafast processes. The predictions of the LLB equations 
related to the linear reversal path for the magnetization 
dynamics [3] as well as to the critical slowing down of the 
relaxation times at high laser pump fluency [5] have been 
confirmed experimentally. 

In ferrimagnetic GdFeCo alloys not only the longitudi- 
nal change of magnetization but also a controllable opti- 
cal magnetization switching has been observed, and this 
has stimulated a great deal of effort to attempt on many 
levels to explain this process, see review in Ref. [TT] . 
The ferrimagnetic materials consist of at least two an- 
tiferromagneticaJly coupled magnetic sublattices. The 
magnetic moments of each sublattice are different, lead- 
ing to a net macroscopic magnetization M(T) defined 
as the sum of magnetization coming from each sublat- 
tice. The main feature of the ferrimagnetic materials is 



that at some temperature, called magnetization compen- 
sation temperature Tm, the macroscopic magnetization 
is zero M(Tm) = 0, although the magnetization of each 
sublattice is not. The angular momentum compensation 
temperature, at which the total angular momentum Ta 
is zero is also of interest. Simplified considerations of the 
ferromagnetic resonance of two-sublattice magnets [T2] 
predict that at this temperature the effective damping is 
infinite and this stimulated investigation of the magneti- 
zation reversal when going through angular momentum 
compensation point [T31 E3] . 

Recently, K. Vahaplar et al. [3|, suggested that the op- 
tically induced ultrafast switching in GdFeCo involves a 
linear reversal mechanism, proposed theoretically in Ref. 
[15] . This is an especially fast mechanism since it is gov- 
erned by the longitudinal relaxation time, which can be 
two orders of magnitude faster than the transverse relax- 
ation time governing processional switching. The mod- 
eling of Ref. [3j was based on macrospin LLB approach, 
essentially treating a ferrimagnet as a ferromagnet. The 
model showed that in order to have the magnetization 
switching a strong field around 20 T was necessary. This 
field can, in principle, come in the experiment with cir- 
cularly polarized light from the inverse Faraday effect. 
More recently, T. Ostler et al. [16] used a multi-spin 
atomistic approach based on the Heisenberg model show- 
ing that the switching occurs without any applied field 
or even with the field up to 40 T applied in the oppo- 
site direction. The predictions for the heat-driven re- 
versal were confirmed in several experiments in magnetic 
thin films and dots using linearly polarized pulses. More- 
over, I. Radu et al. |17j used the same atomistic model 
for the magnetization dynamics to simulate GdFeCo and 
compared the simulation results to the experimental data 
measured by the element-specific x-ray magnetic circular 
dichroism (XMCD). They unexpectedly found that the 
ultrafast magnetization reversal in this material, where 
spins are coupled antiferromagnetically, occurs by way of 
a transient ferromagnetic-like state. 

The latter experiments demonstrate the deficiency in 
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application of the macrospin ferromagnetic LLB model to 
the description of the ultrafast dynamics in a ferrimag- 
netic material GdFeCo. It is clear that the situation of a 
ferromagnetic-like state in a ferrimagnetic material can- 
not be described in terms of a macrospin LLB equation 
in which a ferrimagnet is essentially treated as a ferro- 
magnet. In a ferromagnetic LLB equation the sublattices 
cannot have their own dynamics and thus the processes 
such as the angular momentum transfer between them 
are essentially ignored. In this situation the only pos- 
sible reversal mode is the linear relaxation requiring a 
strong applied magnetic field as was the case of Ref. [1] . 

On a general basis, atomistic models are convenient to 
model ferrimagnetic materials but for modeling of larger 
spatial scales, a macroscopic equation similar to ferro- 
magnetic LLB equation is desirable. This will open a 
possibility to a correct micromagnetic modeling of ferri- 
and antiferromagnetic nano and micro structures at ul- 
trafast timescales and and/or high temperatures. Addi- 
tionally, this can also allow more correct understanding 
of longitudinal relaxation in two-component (for exam- 
ple, ferrimagnetic) compounds, taking into account the 
inter-sublattice exchange. 

In this article we derive a macroscopic equation for 
the magnetization dynamics of a two-component system 
valid at elevated temperatures in the classical case. As 
a concrete example, we consider the disordered GdFeCo 
alloy, the cases of two-component ferromagnets as well as 
ordered ferrimagnets and antiferromagnets can be easily 
deduced. Fig|T]shows a sketch of an atomistic model for a 
ferrimagnetic material and the corresponding micromag- 
netic approximation. The atomistic model is based on 
the classical Heisenberg model for a crystallographically 
amorphous ferrimagnetic alloy |18| and the Langevin dy- 
namics simulations of a set of the Landau-Lifshitz-Gilbcrt 
(LLG) equations for localized atomistic spins. In the 
macroscopic approach each sub-lattice is represented by a 
macrospin with variable length and direction. We use the 
mean field approximation (MFA) to derive a macroscopic 
equation of motion for the magnetization of each sublat- 
tice. It contains both transverse and longitudinal relax- 
ation terms and interpolates between the Landau-Lifshitz 
equation at low temperatures and the Bloch equation at 
high temperatures. We investigate the signs of the relax- 
ation rates of both transition (TM) and rare-earth (RE) 
metals as a function of temperature. We conclude that 
it is a good starting point for performing large scale sim- 
ulations in multi-lattice magnetic systems as the LLB 
equation is for ferromagnetic materials [5J 119) . 



II. ATOMISTIC MODEL FOR A DISORDERED 
FERRIMAGNET. 

The models for binary ferrimagnetic alloys of the type 
A X B\- X , randomly occupied by two different species (A 
and B) of magnetic ions have been previously extensively 
investigated theoretically [20-22 . In such models A and 
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FIG. 1: (Left) Sketch of atomistic regular ferrimagnetic lat- 
tice. Each point represents a magnetic moment associated 
with an atomic site. Magnetic moments of blue points are 
pointing downwards and red ones upwards. (Right) A macro- 
scopic view of partial average magnetization tua = {sa} and 
rriB = (sb) by two macrospins in each sublattice as described 
by the Landau-Lifshitz-Bloch equation. 



B ions have different atomic quantum spin values Sa and 
Sb (Sa ^ Sb)- In the present article we use the classical 
counterpart of these models by considering the classical 
spins with magnetic moments fiA 7^ Ms- We denote A 
specie as TM and B specie as RE. A further but non 
essential simplification is to assume that the interactions 
between spins in the disordered binary alloy are of the 
Heisenberg form with the exchange interactions different 
for different pairs of spins (AA, BB or AB). 

Let us start with the model for a ferrimagnet described 
by the classical Hamiltonian of the type 
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where Af is the total number of spins, (i, j) are lattice 
sites, Hi is the magnetic moment located at lattice site 
i. The external applied field is expressed by H. The 
anisotropy is considered as uniaxial with Di being the 
anisotropy constant of site i. The third sum is over all 
nearest and next-to-nearest neighbor pairs and we have 
considered unit length classical vectors for all lattice sites 
|s.i| = 1. Heisenberg exchange interaction parameter be- 
tween adjacent sites is Jy = Jaa(bb) > if both sites 
(i,j) are occupied by A(B) type magnetic moments and 
Jij = Jab < if the sites (i, j) are occupied by A and 
B respectively. We consider that the ordered TM al- 
loy is represented by the fcc-type lattice. To simulate 
the amorphous character of the TM-RE alloy, x ■ 100% 
lattice sites are substituted randomly with RE magnetic 
moments. 

The magnetization dynamics of this model interact- 
ing with the bath is described by the stochastic Landau- 
Lifshitz- Gilbert (LLG) equation 



Si = ji [si x Hi 



■ Ci\ - 7i^i[Si X [s; X Hj 
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where Xi is the coupling to the heat bath parameter and 
7i is the gyromagnetic ratio. In what follows and for sim- 
plicity we use the same values for TM and RE, 7 TM = 
7 RE = 7 = 1.76 • 10 7 rad s _1 Oc _1 , A TM = A RE = A = 0.1. 
The stochastic thermal fields Ci are uncorrelated in time 
and on different lattice sites. They can be coupled to 
different heat baths (via temperature of phonon or elec- 
tron) and could have different strength of coupling (via 
Xi and Hi) for each atom type (A or B). The correlators 
of different components of thermal field can be written 
as: 



<Ci,a(*)0,/3(*')> = 2 ^ 1 -5 l3 5 al3 5{t - if) (3) 

Hili 

where a, ft are Cartesian components and T is the tem- 
perature of the heat bath to which the spins are coupled. 
The effective fields are given by 
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non- interacting spins pQ. The FPE for the distribution 
function of an ensemble of interacting spins can be de- 
rived in the same way as in the ferromagnetic case [T]. 
The FPE has as the static solution the Boltzmann distri- 
bution function f ({sj}) oc exp [—(3% ({sj})], where W is 
given by Eq. ([!]) and f3 = l/^ksT). Since the exact solu- 
tion is impossible even in the simple ferromagnetic case, 
then, we resort to the mean field approximation (MFA) 
with respect to spin-spin interactions and random aver- 
age with respect to disorder configurations. In the MFA 
the distribution function is multiplicative and we can use 
the same strategy as in the ferromagnetic case [T] , we take 
the distribution function /$ of each lattice site i, which 
satisfy the FPE for a non-interacting spin and perform 
the substitution H =>■ (H^ FA )°° nt , where v =TM or RE 
indicates the sublattices. Thus, we start with the para- 
magnetic LLB equation which was derived in the origi- 
nal article by D. Garanin [I] and is equally valid for the 
present purpose and substitute the external field by the 
MFA one in each sublattice. The corresponding set of 
coupled LLB equations for each sublattice magnetization 
rxij, has the following form: 



The particular values for exchange parameters and the 
anisotropy constants (see Table |l| are chosen in such a 
way that the static properties coincide with experimental 
measurements in GdFeCo 1181. 
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Transition Metal (TM) 


2.217 


8.0725 x 10" 24 


4.5 x 10~ 21 


Rare-Earth (RE) 


7.63 


8.0725 x 10" 24 


1.26 x 10~ 21 


TM-RE 






-1.09 x 10~ 21 



TABLE I: Table with parameters of transition metal (TM) 
and rare-earth (RE) compounds. Anisotropy constant 
Dtm(re) is taken equal for both lattices. Exchange parame- 
ters Jtm(re) /per link are taken in order to give correct Curie 
temperature of pure compounds (x = pure TM or x = 1 
pure RE). Antiferromagnetic exchange parameter J aE -TM IS 
chosen so that the temperature dependence of the TM and RE 
sublattices agrees qualitatively with results of XMCD mea- 
surements of static magnetization [18] . 



III. LLB EQUATION FOR CLASSICAL 
FERRIMAGNET 
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where 
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Here £o,u = |£o B (0 = coth(£) — l/£ is the Langevin 
function, 
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describe parallel and perpendicular relaxation, respec- 
tively, A„jv = ^lvX v / Pijl v is the characteristic diffusion 
relaxation rate or, for the thermo-activation escape prob- 
lem, the Neel attempt frequency. 

Next step is to use in Eqs. Q and ([5) the MFA ex- 
pressions. The MFA treatment for the disordered ferri- 
magnet has been presented in Ref. [T5|. The resulting 
expressions for the fields have the following forms: 



A. Equation derivation 

The idea of the two-component LLB model is presented 
in Fig. [T] Namely, our aim is to evaluate the dynamics of 
the macrosopic classical polarization m = (s) conf , where 
the average is performed over temperature as well as the 
microscopic disorder configurations. 

The dynamics of the mean magnetization can be ob- 
tained through the Fokker-Planck equation (FPE) for 
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where Jo,tm — Q^^tm-tmi ^o : r.e — ^^tm-tmi ^q.re tm — 
qzJ TM . RE , Jo,tm-re = xzJ TM _ RE , z is the number of nearest 
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neighbors between TM moments in the ordered lattice, 
x and q = 1 — x are the RE and TM concentrations. 
The field H^, ff contains the external applied and the 
anisotropy fields acting on the sublattice v =TM,RE. 

The equilibrium magnetization of each sublattice m e iy 
within the MFA approach can be obtained via the self- 
consistent solution of the Curie- Weiss equations 



m RE = B (£ RE ) 



^ RE 



**1tm — B (£tm) 



TM 
TM 



(9) 



The resulting equation Q with expressions Q and |8]) 
constitutes the LLB equation for a ferrimagnet and can 
be already used for numerical modeling at large scale 
since in what follows some approximations will be used. 
The use of these approximations is necessary for under- 
standing the relaxation of a ferrimagnetic system from 
theoretical point of view. We will also get the LLB equa- 
tion in a more explicit and compact form. 

We treat the most general case where the continuous 
approximation in each sub-lattice can be used. Basically, 
in the spirit of the MFA approximation, in each sub- 
lattice we treat the k — mode. In order to handle the 
problem analytically we decompose the magnetization 
vector m„ into two components m v = Yl v + r„, where 
11^ is perpendicular to m K , so that it can be expressed 
as II„ = — [m K x [m K x m„]] /m 2 K , and t v is parallel to 
m K , and it can be expressed as t v — m K (m„ • m K ) jrr? K , 
where k =^ v. 

We can shorten the notation by definition of the fol- 
lowing new variable 0„ K 



m o,i^ up to the first order in H" ff v , under the assumption 
Hg Xl) 3> |H" ff J. Similar to the ferromagnetic case, 
from Eq. ([5| we get (see Appendix A) 



m ,„ 



— xa v + B' u Pii» S o-^ 

m,, mi 



[H'eff^ X m A X m ^ 



m v Jn v 



(12) 



substituting this into Eq. Q and repeating the same 
calculations as in the ferromagnetic case we get the fol- 
lowing equation of motion 
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where B v = B v \J3 J 0tV (m„, m K ) m„J depends on the 

sublattice magnetizations (m„, m K ) and the damping pa- 
rameters are: 
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As a consequence, the MFA exchange field 
in Eqs. ^ and ([8| can be written as the sum of the ex- 
change fields parallel and perpendicular to magnetization 
of the sublattice v. 
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B. Temperature dependence of damping 
parameters 

The temperature dependence of the damping param- 
eters is obtained in the first order in deviations of mag- 
netization from their equilibrium value. Note that in 
Eq. (13) all terms are of the first order in the parame- 
ter H'J S v /H EXjll so that the damping parameters should 
be evaluated in the zero order in this parameter. Conse- 
quently, we can use the following equilibrium expression: 



(15) 



where we have defined a new function J 0j „ (m^rrij,) as 

Jq.m = Jq,v + Jo,uk®ku (m. K ,m„), we remark that Jq, u is 
not a constant but a function of both sublattice magneti- 
zation. The exchange field is, therefore, separated in two 
contributions, a longitudinal one h| x v — ( Jq. m / pL v )m u 

and a transverse one Hg X v = ( Jq mk / fi v )H K . 

In the following we will consider that the transverse 
contribution is small in comparison to longitudinal one, 



i.e. |H 
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- Hg X v . We now expand 



where the sign of the second term does not depend on 
the sign of the interlattice exchange interaction, Jq. VRj . 
The effective damping parameters depend on tempera- 
ture T via temperature-dependent equilibrium magneti- 
zation. The temperature dependence of damping param- 
eters (14), normalized to the intrinsic coupling parame- 
ter, are presented in Fig. [2] for a GdFeCo RE-TM ferri- 
magnet and for various concentrations of RE impurities. 

Let us consider some limiting cases. First we consider 
the simplest case of a completely symmetric antiferro- 
magnet (AFM). In the AFM all the relevant parameters 
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are equal for both lattices, they have the same magnetic 
moments /ii = /12 and the same intra-lattice exchange 
parameters Jo,,,, the inter-lattice exchange parameter is 
also the same Jo,^ K = Jq, kv in contrast to our disor- 
dered ferrimagnet. In this case the equilibrium mag- 
netizations as a function of temperature are the same 
ff^e v{T) = m e K (T) and the effective exchange param- 
eter reduces to Jq^ = Jq v + |J( 



e. the sum of 
the two interactions coming from the intra-lattice and 
inter-lattice exchange. The Neel temperature in the MFA 
reads fceT/v = Jo,i//3 and the damping parameters re- 
cover the ferromagnetic type expression 



i/(afm) 
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37V 
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32V 



(16) 



The use of the critical temperature provides an expres- 
sion in which the damping parameters do not depend 
explicitly on the interlattice exchange, the implicit de- 
pendence comes from the change of the Neel tempera- 
ture as the exchange parameter Jq_ uk varies. There is 
a more simple AFM, with nearest neighbor interactions 
only and one inter-lattice exchange parameter Jq, uk , it 
gives the same result as above and exactly the same as 
for the ferromagnet. 

Next interesting case is when one of the three exchange 
parameters can be neglected. We can consider, for exam- 
ple, a negligible exchange between the rare-earth mag- 
netic moments, it is a good approximation if the impurity 
content is low. Then we can write the effective exchange 
as 



JO, RE 



<A),TM m e,TM ~t~ I Jo,TM-RE 

I Jq.re-tm 



1 m e,TM 



J(),TM (17) 

(18) 



In this case the TM damping parameters can be approx- 
imately expressed with the antiferromagnetic or ferro- 
magnetic (Tjv — > Tq) formula ( 16 ) because in the limit 
x — > the Curie temperature of the disordered ferri- 
magnet is close to ksTc = Jo.tm/3 [TB]. The damp- 
ing parameter for the the RE lattice, however, is dif- 
ferent. It strongly depends on the polarization effect of 
the TM lattice on the RE magnetization. In this case 
close to Tq the polarization effect can be expressed us- 
ing the expansion, B £/3, which for this case reads 

thus, J ,re ~ V(3/3). There- 
fore, we have the following expressions 
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This relation becomes quite important above Tc- We 
observe in Fig. [2] that even for quite large amounts of 
RE of 25% and 50%, the above approximation holds quite 
well. 
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FIG. 2: Damping parameters «y( _)_)(#) (normalized to the 
corresponding intrinsic values) for a pure ferromagnet (FM), 
rare earth (RE) component in a GdFeCo ferrrimagnet and 
a transition metal (TM) in a ferrimagnet as a function of 
reduced temperature •& = T/Tc for three different rare earth 
(RE) concentrations x. The blue solid line represents the 
x — limit which corresponds to a pure ferromagnet (FM). 
(Up) The corresponding curves for a 25% concentration of 
RE. (Middle) The corresponding damping parameters for a 
50% alloy. (Bottom) Damping values for 75% RE amount. 
It can be also seen as a RE doped with a 25% of transition 
metal (TM). 



If the inter-lattice exchange is large in comparison to 
the intra-lattice one then the equilibrium magnetization 
of both lattices is similar and the damping parameters 
behave similar to those of the FM damping parameters, 
presented above. This case is in agreement with a con- 
centration of 75% of RE in Fig. [2] (down). As predicted, 
we observe that the damping parameters are very similar 



G 



for both sub-lattices. 

Note that these damping parameters should be dis- 
tinguished from those of the normal modes (FMR and 
exchange) with more complicated expressions which can 
be obtained via linearization of the set of two-coupled 
LLB equations [28], similar to the LLG approach. 



C. Longitudinal relaxation parameters 



The function 1 — B v jm v in Eq. (13 1 is a small quan- 
tity proportional to the deviation from the equilibrium in 
both sublattices. It can be further simplified as a func- 
tion of the equilibrium parameters after some algebra. 
Similar to the ferromagnetic case, the ferrimagnetic LLB 
equation can be put in a compact form using the notion 
of the longitudinal susceptibility. 

The initial longitudinal susceptibility can be evalu- 
ated on the basis of the Curie- Weiss equations ([9]). Let 
us assume that in the absence of an external field, the 
equilibrium sublattice magnetizations m TM and m RE are, 
respectively parallel and antiparallel to the z-axis (a 
stronger condition of the smallness of the perpendicular 
components can be also applied). The z-axis is chosen 
such that it is the easy axis of the magnetic crystal. To 
evaluate the longitudinal susceptibility, the field should 
be applied parallel to the easy direction, then in the ap- 
proximation of small perpendicular components (large 
longitudinal exchange field) we can neglect in the first 
approximation the possible change of directions of m RE 



and m TM . In order to calculate the susceptibility, we 
expand the right-hand side of Eq. ^ in terms of the 
external field: 



iu(T, H z ) « m v {T, 0) + ti„H z pB' v 1 
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where B v = B v (j3(j, v H EX<v ) and its derivative B' v — 
B^PfJiuH-Ex u) are evaluated in absence of applied and 
anisotropy fields. Then, 
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Thus, the longitudinal susceptibility of one sublattice is 
expressed in terms of another: 
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-x K ,\\ 



(23) 



Finally, we obtain two coupled equations for x RE n and 
Xtm,||i solving them, wc get the MFA expression for the 
susceptibilities: 
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The longitudinal susceptibility Xv.\\ therefore, a func- 
tion of temperature which we have called G V (T). It 
tends to zero at low temperature and diverges approach- 
ing Curie temperature Tq of the magnetic system, sim- 
ilar to the ferromagnetic case. The function G v = 
(\Jo,vk\/(1'i>)x'v,\\ can be seen as a reduced longitudinal 
susceptibility. 

Now we derive an approximate expression for the small 
quantity 1 — B v /m v as a function of equilibrium quan- 
tities and the deviation of each sublattice magnetization 
from its equilibrium. In the first approximation, we ex- 
pand the function B v jm v near the equilibrium, as was 
done for the ferromagnet. The function B v in the zero 
order in perpendicular field components, H"** v /Hex,v> 
can be written as a function of m v and m K as follows 



B v sa B v (0[J o , v m v + | Jo, vr \t k ]) (25) 
where t k = \(ra. v -Ta K )\/ m v is the length of the projection 



of the magnetization of the sublattice k onto the sublat- 
tice v. We expand the function B v lm v in the variables 
m„ and m K near the equilibrium : 



rn,. 



B, 



1 dB v 

m v dr K 



1 



( dB u 



:B V 



= l-[l-p\ v B' v 



(26) 



cq 



Joq 



Wo,uk\K] 

here 5m y = m v - m e>v , with m e , v = B v (Pfi u H EXiU ), 
where H EX L , is evaluated at the equilibrium, and St k = 
t k - T e!K , where r GiK = |(m e , y • m e , K )\/m e , u and it corre- 
sponds to the projection of the equilibrium magnetization 
vn. e K onto the other sublattice magnetization direction. 
It is easy to show that dT K /dm v = 0. Similar to the fer- 
romagnetic case, we would like to arrive to a simplified 
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expression as a function of sublattice susceptibilities. For 
this purpose, we divide the above expression by fi u /3B v 



1 - B u /m v 



1 5m v 

Xu,\\ m e , v 

1 Sm„ 



X»,\\ 



Xn,\\ 



(27) 



where we have used Eq. (23) and the function G K 



| Jo,vk\Xk.\\/( j 'i> nas now more sense. Thus, the contribu- 
tion to the dynamical equation Q of the exchange in- 
teraction (the LLB equation with longitudinal relaxation 
only) given by Eq. (27) reads 



Form 1 

Finally, we collect all the above derived approximate 
expressions and we finish up with the compact form of 
the LLB equation for the reduced magnetization vector, 
m v =M v /M v (T = 0K) 

r XT 1 V ' Hoff,!/) 

m v = 7 !/ [m !/ x H cff:Jy J - j„a» = — — m„ 



v [m„ x [m„ x K eS v ]} 



m- 



(33) 



where the effective field H c ff „ for sublattice v is defined 
as 



1v 



\EX 



< ~' K Sm v ~ St k I xa u (28) 



Xv,\\ 



Note that the first term defines the intralattice relaxation 
of the sub-lattice (for example, TM) to its own direction. 
The second term describes the angular momenta trans- 
fer between sublattices driven by the temperature. This 
equation has the form 



(29) 



and it gives the exact LLB equation for the case when 
the average magnetization of the two sublattices remain 
always antiparallel. 



H 



eff,i 



H + H 



IT 



1 

2a: 



m e.i/ / 2A„ K v , ; 



m^34) 

and the relaxation parameters aY, and are given by 



Eqs. (14) 



Or in a more explicit form, as a function of sub-lattice 
magnetizations m„ and its values at the equilibrium 
m„ „: 



m„ = y„[m„ x H eff)I/ ] - ^ v a\ 



m » ■ H cff.„ 



v [m„ x [m„ x Heff, v ] 

7^_L 2 



(35) 



D. Final forms of the LLB equation 

In order to be consistent with the ferromagnetic LLB 
equation (and the Landau theory of phase transitions), 
we expand the deviations 8m v (St k ) around m\ v {t^ v ) 
up to the quadratic terms. Similar to FM case we write: 



<5m„ 



m„ 



( 2 2 

m„ - m„ 



(30) 



Therefore we can write the effective longitudinal fields as 



H 



off. I I 



2A„ 



2A„ K 



where in order to shorten the notations we have defined 
the longitudinal rates as: 



where we have defined the longitudinal field, HjJ ff u , as 



H 



eff.i/ 




and the effective field, H G g v , reads 



H e ff,i/ — H + H A .„ H — ^^m K 



m„ (36) 



In Eq. ( 35 ) also the temperature dependent damping 



(31) parameters are given by Eqs. (14) 



Form 2 



A: 



1 

Xv,\\ 



(1 + G K ), A-i = ^^iwith^ N 



flu 



where G K is also expressed in terms of the longitudinal 



susceptibility via Eq. ( 24 ) 



It is also interesting to put the LLB equation in a more 
symmetric form in terms of the macroscopic magnetiza- 
tion, M„ = x v jjL v xa v /v v , where x v stands for the concen- 
tration of sites of type v =TM or RE (x v — x for RE and 
x v = q for TM), /i„ is the atomic magnetic moment of 



the lattice v and v v is the atomic volume. We multiply Next, we substitute Eq. (23) into Eq. (41 1 and we get the 



each sublattice LLB equation ( 35 1 by the corresponding 
factor, for example, in the case of TM by qfJ, T M/ v TM an d 
we obtain 



M„ = 7 „[M V x H eff)l/ ] - % 
[M„ x [M„ x H cff , 



- L 



M, 2 



(37) 



where the effective hclds read: 



H 



2A„ 



M 2 











' M„- 






M e , K 



M„, (38) 



The rate parameters are A„ K = v v t^ UK j [i v x v and the ef- 
fective field, H c ff jly , has the following form: 



H, 



eff.i 



H + H A , y +ylM K . 



Here the exchange parameter is introduced as A — 
z^tmre/Mre^tm- The damping coefficients Lm,, and 
£j_,„ read 



IV. RELAXATION OF MAGNETIC 
SUBLATTICES 

The rate of the longitudinal relaxation is temperature 
dependent through the parameters such as the damping 
parameters aj|, see Eq. (14) and Fig. |5J and the longi- 
tudinal susceptibilities. The sign of the rate, T v ^ 0, de- 
pends on the instantaneous magnetization values. From 
Eq. ( 28 1 we can consider the following lines separating 



different relaxation signs: 

0m v = — — —rOT K = Xvk,\\ot k , (39) 

where we have defined the dimensionless variable Xvk.\\i 
which describes the effect of the change in one sublat- 
tice on the other. This variable can be interpreted as a 
susceptibility Xvk\\ = 5ra v /5m K . Indeed, we can expand: 

m v (T,5m v ,5m K ) « m v {T, 0, 0) + pjQ <v B' v 8m v 

+/3\J ,„ K \B'Jm K (40) 

Now using that by definition 5m u = m v (T, 5m v , 5m K ) — 
m v (T, 0, 0) , we obtain 

^ = ^{t^wJ (41) 



relation between the susceptibilities, exactly described by 
Eq. p9l. 



The problem of relaxation sign is, therefore, reduced to 
the study of the sign of the function 8m v — Xvk,\\& t k- Let 
us assume the equilibrium state that is close to Tc, de- 
scribing the situation during the ultrafast laser-induced 
demagnetization [T7]. Fig. |3] shows three possible instan- 
taneous rates for T = 0.95Tc, depending on the relative 
state of both sublattice magnetizations. The lines sep- 
arating different relaxation types are straight lines with 
the slope Xu K ,\\{T). 

In the following we use atomistic LLG Langevin sim- 
ulations described in Sec. [TT] as well as the integra- 
tion of the LLB equation Q for the same material pa- 
rameters, see Table [TJ In order to compare MFA based 
LLB equation and the atomistic simulations, we have 
re-normalized exchange parameters, as described in Ref. 
[181 . In the atomistic simulation the system size is taken 



as N = 60 3 



i.e. 



37V coupled differential equations has to 



be solved simultaneously within this approach, whereas 
only 6 (two sublattices and three components for each) in 
the macrospin LLB approach. We compare the different 
relaxation regions depending on the instantaneous mag- 
netic state with those predicted by the LLB equation and 
depicted in Fig. [3| The initial conditions in the simula- 
tions are the following: in all three cases we start from 
a an equilibrium state at T = 600 K (for the considered 
concentration x — 0.25 we get Tc — 800 K). After that 
for the situations of FigQa) and (e) we put one of the 
sublattice magnetizations equal to zero, m TM(RE) = 0. In 
the atomistic approach this is done by totally disordering 
the system. Finally, the temperature is set to T — 0.95Tc 
and the relaxation of both sublattices is visualized. The 
results are presented in Fig. |4j 

For the region m RE 3> m TM above the green line in 
Figji] the rate for the TM is positive, T TM > 0, thus 
the TM magnetization will increase while r RE < and 
the RE magnetization will decrease. Thus, we have ini- 
tially a dynamical polarization of TM by RE. As it can 
be seen in FigQa),(b) initially the TM magnetic order 
increases from a totally disordered state, while the RE 
relaxes directly to the equilibrium, i.e. the sign of the 
RE rate is always the same. In the central region of 
Fig(3j between green and red lines, both magnetizations 
go to the equilibrium by decreasing their value, sec Fig. 
|4^c)(d). Finally, in the low region of Fig [3] the situation 
is symmetric to the upper region but now TM magne- 
tization decreases and the RE magnetization increases 
initially, see Fig. |4|e)(f). Thus, the predictions of the 
LLB equation are in agreement with full atomistic simu- 
lations which also provides a validation for our analytic 
derivation. 

As a representative example, in GdFeCo near the mag- 
netization reversal the situation is the following [17] : the 
TM magnetization is almost zero, m TM ~ and the 
RE has finite magnetization value m TM > 0. This hap- 
pens due to the fact that the Gd sublattice is intrin- 
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Ctm < 0, 

r RE > o 



0.1 0.2 0.3 0.4 0.5 0.6 

m T M 



FIG. 3: Different longitudinal relaxation regions for T/Tc 
0.95 for parameters of the GdFeCo alloy with x = 0.25. 



sically slower than the FeCo one due to a larger mag- 
netic moment. This situation corresponds to the up- 
per region in Fig. [3] where the rates are r TM > and 
r RE < 0. Under these circumstances the RE magne- 
tization dynamically polarizes the TM sublattice mag- 
netization through the interlattice exchange interaction 
#ex,tm-re ~ |^o,tm-re|"i re > 0. Consequently, the TM 
magnetization goes opposite to its equilibrium position 
m™ = [see Fig. |4ja)-(b)]. The existence of opposite re- 
laxation signs in TM and RE is consistent with a recently 
reported ferromagnetic state in a ferrimagnetic materials 
Ref. [T7], however it does not necessary lead to it. Nor 
it necessary means the switching of the TM magnetiza- 
tion, as was suggested in Ref.pl]. To have a switching 
one should cross the line m™ = which cannot be done 
within the approach of longitudinal relaxation only which 
only describes the relaxation to the equilibrium. The 
crossing of the line m™ = can be only provided by a 
stochastic kick which is always present in the modeling 
using stochastic atomistic approach [I|)J |T7] . This topic 
will be the subject of future work. 



V. THE LLB EQUATION AND THE 
BARYAKHTAR EQUATION 

In this section we would like to discuss the differences 
between the LLB equation and the equation derived by 
V. Baryakhtar 25j and used in Ref. 24 to explain 
the ultrafast magnetization reversal and the transient 
ferromagnetic-like state in ferrimagnets. The Baryakhtar 
equation was derived from the Onsager principle which 
in general is valid near the thermodynamic equilibrium 
only. The general derivation is based on the symme- 
try approach. Another strong supposition made in its 
derivation is the separation of the timescales: the ex- 
change interaction timescalc and the relativistic interac- 
tion timescale (defined in our case by the parameter A) 
are assumed to be separated. The resulting equation has 
the following form: 



LLG-Langevin 



LLB 




(b) 
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0.0 
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-0,2 
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t 
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(c) 
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(e) 
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0.1 






0.0 1 






r 




-0.1 







(f) 
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time [ps] 



FIG. 4: Comparison between atomistic LLG-Langevin and 
macrospin LLB calculations of the longitudinal relaxation of 
the GdFeCo alloy (x — 0.25) corresponding to the three differ- 
ent relaxation cases in Fig[3] In the left column we show atom- 
istic LLG-Langevin multispin simulations and in the right 
one- the LLB macrospin calculations. The graphs (a) and 
(b) correspond to the region with T T m > and r RE < 0. The 
graphs (c) and (d) correspond to the region with F T m < and 
r RE < 0. The graphs (e) and (f ) correspond to the region with 
f T M < and f RE > 0. 



1 dM v 
7„ dt 



A e (H v — H K ) + \ V H V 



(42) 



Here v — TM, RE, A„ describes transfer of the angu- 
lar momentum from sublattices to the environment, A e 
is of the exchange origin and stems from spin-spin in- 
teractions, conserving the total angular momentum but 
allowing for the transfer of angular momentum between 
the sublattices. The effective fields defined as H,, = 
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— SW/SM y are derived from the magnetic energy W. In 
Ref. [24] the authors used the Landau type free energy 
expansion near the critical temperature, corresponding 
to the form Eq. p0| . 

In comparison to the Baryakhtar equation, the LLB 
equation, derived here includes the transverse exchange 
mode and allows the transfer of the energy or momentum 
between the longitudinal and transverse motion. The fcr- 
rimagnetic LLB equation has three terms among which it 
is the precession term which conserves the total angular 
momentum. The precession in the interlattice exchange 
field given by [m TM x m RE ] allows the transfer of an- 
gular momentum between sublattices. The longitudinal 
and transverse relaxation terms which are related to the 
coupling to the heat bath are both proportional to A. 
Differently to ferromagnets, both the transverse motion 
given by precession and transverse relaxation terms are 
not negligible on the femtosecond timescale in compari- 
son to longitudinal motion because in both cases the field 
acting on both motions is of the exchange origin. 

In principle the ferrimagnetic LLB equation can be cast 
in a form, similar to the Baryakhtar equation if we re- 
strict ourselves to longitudinal motion only, considering 
the antiparallel sublattices alignment. For the longitudi- 
nal relaxation only (see Eq. (28)) we have the following 
expression 



(43) 



where H' v = —(^ bL )m'^/m l ,, stands for the fields coming 
from interaction of each lattice with itself and H' K -with 
the opposite sublattice. One can see that the sign of the 
effective field coming from the other sublattice is opposite 
for the LLB Eq. ( 43 1 and the Baryakhtar equation Eq. 
In order to illustrate the consequence of this, we 



(421. 



can compare the equations for the limiting case close to 
Tc- In this case the Baryakhtar equation (see Eq. (1.33) 
in Ref. reads: 



(44) 



where m u z is the absolute value of the z-component of 
the magnetization in the sub-lattice v and we explicitly 
considered that the sign of ^-components is opposite for 
the sublattice v and k. In the same limit, considering 
to T m(re) = TOe,TM(R E) + <5m TM(aE) , and following Eq. ( [28] ) 
the LLB equation takes a similar form: 



Xu,\\ 



Xu,\\ 



(45) 



Note that for the LLB equation the contribution of the 
opposite sublattice is negative while for the Baryakhtar 
equation it is positive. This has important consequences 



in the longitudinal inter-lattice relaxation of the sub- 
lattices, changing the results of Fig. [3] 

In Fig. [5] we show the temperature dependence of the 
ratio of partial susceptibilities, Xk,\\/Xv,\\ appearing in 
Eq.(|4"5"l). We can see that at temperatures not very close 
to Tc'- Xtm,||/Xrb,|| 1 and the contrary behavior close 
to Tc- Thus for the TM and temperatures close to Tq the 



second term in the r.h.s. of Eq.(45) could be neglected 
and the third term with the opposite sign can compete 
with the first one, leading either to slowing down of the 
relaxation rate or even to changing its sign, as presented 
in Fig. |4|a) and (b). The behavior of RE on the contrar- 
ily is dominated by this term and the sign of relaxation 
cannot be changed, as is seen in the same figure. Obvi- 
ously this behavior cannot be described by Eq. ( 44 ) where 



all terms have the same sign. In order to have the oppo- 
site relaxation sign, one has to assume for this equation 
a priori that the signs of the z-components of magnetiza- 
tion in both sub-lattices are the same, i.e. to start with 
the ferromagnetic-like state without specifying its origin. 
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FIG. 5: Temperature dependence of the ratio between longi- 
tudinal susceptibilities for parameters of the GdFeCo alloy. 



Finally, we would like to note that because we have 
treated the spin-spin interaction in MFA we have lost 
correlation contribution. Consequently, both LLB and 
Baryakhtar equations do not describe the energy trans- 
fer from the uniform modes into nonlinear spin waves and 
vice versa. In ferromagnets |26j this contribution is usu- 
ally two or three orders of magnitude smaller than the 
contribution to relaxation through the coupling to the 
bath. At this stage we do not know how large this contri- 
bution can be in ferrimagnets. In Ref. |26j the contribu- 
tion of nonlinear spin waves was artificially incremented 
by using a random anisotropy to cause non-coliniarities. 
In principle, in ferrimagnets one can see a small amount 
of RE as precursor of non-coliniarities, with the strength 
of the order of interlattice exchange parameter J TM re- 
For completeness, a microscopic treatment of the spin 
wave contribution would be desirable, we let this task for 
the future. 
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VI. CONCLUSIONS 

We have derived the Landau-Lifshitz-Bloch equation 
for a two-sublattice system such as a GdFeCo ferrimag- 
net for which an ultrafast switching has been reported 
[T41, [T7] . Although in our derivation we refer to a TM- 
RE alloy, it is equally valid for a two-component ferro- 
magnet, as well as for an antiferromagnet. The general- 
ization to more components is straightforward. The new 
equation constitutes an important step forward in clas- 
sical description of the dynamics of ferrimagnets which 
is traditionally based on two-coupled macroscopic LLG 
equations. For example, the FMR and exchange modes 
have recently attracted attention due to possibility to 
optically excite them [T31 H7j. Their temperature de- 
pendence can be now correctly understood in terms of 
our approach |28j . Furthermore, recent ultrafast dynam- 
ics experiments using XMCD showed different sublattice 
dynamics on ultrafast timescale in a two-sublattice mag- 
nets such as GdFeCo [UJ or FeNi [29] , which can be mod- 
eled using this new approach. Finally, this equation can 
serve in the future as a basis for multiscale modeling in 
two-component systems at high temperatures and/or ul- 
trafast timescales, the same way as the LLB equation 
for ferromagnets [T5]. This also opens a possibility for 
micromagnetic modeling of ultrafast dynamics in large 
structures, such as sub-micron and micron-size ferrimag- 
netic dots, whose dimensions do not allow modeling by 
atomistic approach. Similarly, it will be useful for static 
micromagnetic modeling at high temperatures, such as 
thermally-driven domain wall motion in nanostructures. 

The LLB equation correctly shows the possibility to 
reverse the sign of relaxation at high temperatures and, 
therefore, is consistent with the existence of a recently 
reported ferromagnetic state in a ferrimagnct 17\. The 
validity of the approach has been checked against full- 
scale atomistic simulations presented in Fig. |4j However, 
unlike the equation, derived by Baraykhtar and used re- 
cently to describe the GdFeCo switching [24], it is not 
based on the separation of timescales and on the On- 
sager principle. Instead, both the coupling to the exter- 
nal bath and the exchange interaction form part of the 
same longitudinal and transverse relaxation terms. We 
show important differences in the resulting form of the 
equation. 

Unfortunately, at the present time the compact deriva- 
tion was possible only under some assumptions. The 
employed conditions certainly allow to describe the nor- 
mal modes such as ferromagnetic resonance and anti- 
ferromagnetic exchange precessional modes in ferrimag- 
nets |28j . The same way the approximation is suffi- 
cient to describe the switching of ferrimagnet if it occurs 
through a linear reversal path [U [M] or if sublattices 
non-collinearities are not too large. Weather the applied 
approximation completely describes the situation of the 
ultrafast reversal is an open question which we will in- 
vestigate in the future. For modeling, the initial param- 
agnetic equation Q with the MFA field Q and pi can 



always be used, providing the check for the approxima- 
tion. Finally, up to now we were not able to derive a 
compact expression for the equation above Tq which is 
also a necessary step for the full modeling of the ultrafast 
switching. 
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Appendix A 

In this appendix we present detailed derivation of Eq. 
( [12] ). We start from Eq.([5|: 



m ,, = Bfa iU )H V) £ 0jt/ = (H^ FA ) co "\ (Al) 



H cff,i/- 



where u„ = €o,*/&,* and (H*"*)-* = H EXji 
Here H" ff contains the anisotropy, applied and the per- 
pendicular component of the exchange field (see section 
III. A). In the case of a strong homogeneous exchange field 



H 



EX., 



> 



H cff,i 



first order in H' cff v 



the MFA field can be expanded up to 
as 



Till . H" 



EX.i/ 



H 



EX,i/ 



(A2) 



MFA \conf 



(Hf A ) 



can be writ- 



H ex,i/ • H eff,J / H lx,v Expanding the Langevin 



Therefore, £ ,i/ = 
ten as £ ,i/ = Cex.i/ + ^ witn £ex,i/ > 
where we identify £ex,i/ = P^fB EX v an< ^ = 

PHu ( 

function around £ex,i/ we get 

B (£„,„) -B v + B'Ji 

and 

/ 



(A3) 



H 



EX., 



H' 



EX,v 



V 



h h" 

"■EX.u n cff.i 



( H EX,J) 



,(A4) 



wher e B v = B (£ex, ^) a n d B' v = B' (C ex.i/)- Substituting 
Eq. (lA3| and Eq. (|A4| in Eq. ([All and neglecting the 



terms quadratic in H" s ^/I^ex v \ we § et 



mo,. 



B„ 



H 



ex,, 



H 



eff.i 



n EX.i/ ' n eS,u ) EX,i/ 



EX., 



EX., 



(h| x ,i/ ' Hg ff l/ ) H EX 



( H EX,v) 



(A5) 
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Using the vector calculus identity (ax b) x c = b (a • c) 
a (b • c) Eq.( A5 1 can be written as 



nifl.i/ — B v __f X '" + B' u P(j, u 



H 



EX„ 



EX,i/ 



B v 



n eflE,y x EX.i/ 



x H 



EX, 



EX.i 



(4x,„y 



(A6) 



Finally, we use h| x >u = (j^u/fJ-v) [see Eq. (JTTJ)] 
Eq. (|A6| and obtain Eq. (fl2l 



B v Hv 
m v J , v 



r" 



m » ' H cff^ ) m 

ml 

x m, 
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